ON THE PERIODIC LORENTZ GAS 
AND THE LORENTZ KINETIC EQUATION 

FRANCOIS GOLSE 

Abstract. We prove that the Boltzmann-Grad limit of the Lo- 
rentz gas with periodic distribution of scatterers cannot be de- 
scribed with a linear Boltzmann equation. This is at variance 
with the case of a Poisson distribution of scatterers, for which 
the convergence to the linear Boltzmann equation was proved by 
Gallavotti [Phys. Rev. (2) 185, 308 (1969)]. The arguments 
presented here complete the analysis in [Golse-Wennberg, M2AN 
Model. Math, et Anal. Numer. 34,1151 (2000)], where the impos- 
sibility of a kinetic description was established only in the case of 
absorbing obstacles. The proof is based on estimates on the distri- 
bution of free-path lengths established in [Golse-Wennberg loc.cit.] 
and in [Bourgain-Golse-Wennberg, Commun. Math. Phys. 190, 
491 (1998)], and on a classical result on the spectrum of the linear 
Boltzmann equation which can be found in [Ukai-Point-Ghidouche, 
J. Math. Pures Appl. (9) 57, 203 (1978)]. 
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1. Introduction 



About 100 years ago, Lorentz [15] proposed the following linear ki- 
netic equation to describe the motion of electrons in a metal: 

(1.1) (d t + v ■ V x + ±F(t, x) ■ V v )f(t, x, v) = N at r 2 at \v\C(f(t, x, •))(«) 

where /(£, x, v) is the (phase space) density of electrons which, at time 
t, are located at x and have velocity v. In Eq. (11.11) . F is the electric 
force field, m the mass of the electron, while N at and r at designate 
respectively the number of metallic atoms per unit volume and the 
radius of each such atom. Finally C(f) is the collision integral: it acts 
on the velocity variable only, and is given, for all continuous <p = < P{ V ) 
by the formula 

(1.2) C((j))(v)= (<p(v-2(v-w)u)-<f>(v))caa(v,u)du. 

J |u>|=l,1Mt>>0 
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In the case where F = 0, Gallavotti [HI [12] proved that Eq. (11.11) de- 
scribes the Boltzmann-Grad limit of a gas of point particles undergoing 
elastic collisions on a random (Poisson) configuration of spherical ob- 
stacles. His result was successively strengthened by Spohn [18], and by 
Boldrighini-Bunimovich-Sinai [I]. 

In the presence of an external, non-zero electric force F and for 
the same random configuration of absorbing obstacles as in [TTl [12], 
Desvillettes-Ricci JTU] proved recently that the Boltzmann-Grad limit 
of a gas of point particles leads to a non-Markovian equation — see 
also an earlier, similar observation by Bobylev-Hansen-Piasecki-Hauge 
i- 

The case of periodic configuration of obstacles, perhaps closer to 
Lorentz's original ideas, completely differs from the random case. In 
the case of absorbing obstacles, and without external force F, several 
results suggest that the Boltzmann-Grad limit is non-Markovian [8l[T4"t 
[9]. However, all these results are based on explicit computations that 
are possible only in the case of absorbing obstacles. 

In the present note, we show that neither equation (11 .ip nor any 
variant thereof can describe the Boltzmann-Grad limit of the periodic 
Lorentz gas with no external force (F = 0) and in the case of reflecting 
obstacles. 

The crucial observation (already made in [8j[Tl]) is that the distri- 
bution of first hitting times satisfies, in the periodic case, an inequality 
stated below as Theorem 12.11 that prevents these first hitting times 
from begin exponentially distributed. The probabilistic representation 
of equation (II. ip by a jump process with exponentially distributed jump 
times plus a drift (see for instance [16J) suggests that the Boltzmann- 
Grad limit of the periodic Lorentz gas with reflecting obstacles cannot 
be described by (11.11) . 

In the present paper, we give a complete proof of this fact (see The- 
orem UJ] below), which by the way does not appeal to the probabilistic 
representation of (II. ip . In addition, the method of proof used here pro- 
vides an explicit estimate of the difference between the single particle 
phase space density for the periodic Lorentz gas in the Boltzmann-Grad 
limit and the solution of the Lorentz kinetic equation (II. ip — see the 
inequality (14.181) and the discussion thereafter. This result has been 
announced in [T3] 

Recently, Ricci and Wennberg [17] have considered the following in- 
teresting variant of the periodic Lorentz gas studied here. Their mi- 
croscopic model consists of gas of point particles in a periodic con- 
figuration of obstacles that are bigger than in the Boltzmann-Grad 
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scaling considered in the present paper. However, some of these ob- 
stacles are removed with a probability carefully chosen in terms of the 
obstacle radius so that the mean collision time remains of order one. 
They proved that the expected single-particle phase-space density so 
obtained converges to a solution of the Lorentz kinetic equation (11. ip . 
This result suggests that the failure of (11.11) to capture the Boltzmann- 
Grad limit of the periodic Lorentz gas is a very unstable phenomenon 
that is specific to the periodic setting and likely to disappear whenever 
some amount of randomness is injected in the microscopic system. 

2. The periodic Lorentz gas 
Let D G N, D > 2. For each r G (0, |), we consider the domain 
Z r = {x G R D | dist(x, Z D ) > r} , 

which is usually referred to as "the billiard table" . 

The Lorentz gas is the dynamical system corresponding to a cloud 
of point particles that move freely in Z r — collisions between parti- 
cles being neglected — and are specularly reflected at the boundary 
of each obstacle — the obstacles being the balls of radius r centered 
at the lattice points, i.e. the connected components of Z°. Since col- 
lisions between particles are neglected, one can equivalently consider 
a single particle whose initial position and velocity are appropriately 
distributed in Y r x S D_1 . 

In this dynamical system, a particularly important notion is that 
of "free path length": see the review article by Bunimovich |5j, pp. 
221-222. 

The free path length — or "(forward) exit time" - - for a particle 
starting from x G Z r in the direction v G S D_1 is defined as 

(2.1) r r (x, v) = inf{t > | x + tv G dZ r } . 

The function r r is then extended by continuity to the non-characteristic 
part of the boundary of the phase-space, i.e. to 

{(x, v) G dZ r x S D_1 \ v -n x ^0}, 

where n x designates the inward unit normal field on dZ r (pointing 
toward Z T) i.e. outside of the scatterers). Because Z r is invariant 
under Z D -translations, one has 

r r (x + k, v ) = r r (x, v ) for each (x, v) G Z r x S D_1 and k G Z D . 

Hence r r can be seen as a [0, +oo]-valued function defined on Y r x S D_1 
(and a.e. on Y r x S 1 ^ 1 ), where Y r = Z r /Z D . 
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Figure 1 . The billiard table Z r and the punctured torus 

Yr 



Whenever the components of v e S D_1 are rationally independent 
- i.e. if k • v ^ for each k 6 Z D \ {0} — each orbit of the linear 
flow x i— > x + tv is dense on T D = H D /Z D , and thus r r (x,v) < +oo 
for each x G Z r . 

On the measurable space Y r x S D_1 equipped with its Borel a-algebra, 
we define /i r as the probability measure proportional to the Lebesgue 
measure on Y r x S ' 1 : in other wor d£] 

dydv 



dfir(y,v) 



\YJ\S 



D-ll 



1 If A is a measurable d-dimensional set in R D (d < D), we denote by |A| its d- 
dimensional volume. Here \Y r \ is the Lebesgue measure of any fundamental domain 
of the quotient space Y r , i.e. of the unit cube with a ball of radius r removed. 
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Define the distribution of r r under \i r by the usual formula 
$ r (t) := ^ ({(y, v)eY r x S^ 1 | T r (y, v) > t}) . 
Theorem 2.1. Let D > 2. There exist two positive constants C\ and 



D-1 



C*2 such that, for each r G (0, \) and each t > 1/r 

In the theorem above, the lower bound in the case D = 2 and the 
upper bound for all D > 2 were proved by Bourgain-Golse-Wennberg 
[8|; the lower bound was extended to the case of any D > 2 by Golse- 
Wennberg [14j . More precise results concerning $ r (t/r) in space di- 
mension D = 2 have recently been obtained by Caglioti-Golse [9] and 
Boca-Zaharescu [3]; however, only the result above (Theorem 12.11) is 
used in the present paper. 

3. The linear Boltzmann equation 

The linear Boltzmann equation for a fred gas of particles moving at 
speed 1 is 

(3. 1) (d t + v V x )f(t, x, v) = C(f(t, x, -))(v) 

where / = f(t,x,v) is the single particle phase-space density, also 
known as "distribution function". In other words, f(t,x,v) is the den- 
sity of particles which, at time t, are located at x e R D and move in the 
direction v G S D_1 . The term C(f) is the collision integral; specifically, 
C is an integral operator of the form 



(3.2) C((j))(v)=a k(v,w)((j)(w) - (j){v))dw 

Js - 1 

where a > and k G C(S D_1 x S 15-1 ) satisfies the following properties 
(3.3) 

k(v , w) = k(w, v ) > , / k(v , w)dv = 1 , for all v , w G S ^ 1 . 



Changing variables according toa;i— > w = v — 2(v ■ uo)uj in the Lorentz 
collision integral for D = 3 shows that 

v — 2(v ■ uj)uj) — 4>(v)) cos(f , uj)duj 

k(v,w)(4>(w) — 4>{y))dw with k{y,w) = ^- . 

s 2 



2 I.e. without external force. 
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Hence, the Lorentz kinetic model fll.ip - fll.2B is a special case of the 
linear Boltzmann equation fl3.ip - fl3.2ft . 

In the sequel, we shall restrict our attention to the case where the 
linear Boltzmann equation is posed in a periodic box. Without loss of 
generality, we assume that this periodic box has size 1, and that S D_1 
is endowed with its rotationally invariant unit measure, henceforth de- 
noted by dv. Finally we denote by (•) the average with respect to both 
variables x and v: 



(3.4) (0) = / / cj)(x } v)dxdv. 

Consider then the unbounded operator A on L 2 (T D x S D_1 ) defined 

by 



(A(J))(x, v) = — v ■ V x (f)(x, v ) — a(j)(x, v ) + a / k(v,w)4>(x,w)dw 

Js - 1 

with domain 

D(A) = {(f) E L 2 (T D x S^ 1 ) | v ■ V x (f) G L 2 (T D x S^ 1 )} . 

We recall the following result, originally proved in [TH] for the more 
complicated case of the linearization of Boltzmann's equation at a uni- 
form Maxwellian state. 

Theorem 3.1. Let a > and k G C(S D_1 x S D_1 ) satisfy the as- 
sumptions A3.3\) . Then the operator A generates a strongly continuous 
contraction semigroup on L 2 (R D x S ^ 1 ), and there exists positive con- 
stants c and 7 such that 



-tA. 



e <p- W||i,2(to x so-i) S ce ' ||</>||l2 (T £> xS £>-i) 



< Ce ^WdW r2^C vS C-n , t>0 



for each G L 2 {T D x S^ 1 ) 

4. The non-convergence result 

Throughout this section, we denote by e the sequence (l/n) ne rsr*. For 
some given r* G (0, |), we set r = r^e 1 ^ 13-1 ), and we denote by f2 e the 
open subset eZ r with this particular choice of r. 

4.1. Reflecting vs. absorbing obstacles. First we define the bil- 
liard flow on the scaled billard table fl e . It is a one-parameter group 
on Q e x S D_1 denoted by 

(4.1) S e t : Vl € x S - 1 3 (x,v) i-> (X^(x,v),V t e (x,v)) G fi e x S^ 1 
and defined in the following manner: 
• Sq(x,v) = (x,v), 
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Figure 2. The billiard dynamics on fi e 



• if t is such that X^(x,v) ^ dfl t , then S^(x, v) is smooth in the 
time variable near that particular value of t and one has 

(4.2) jXt{x, v) = V7(x, v) , jVf(x, v) = ; 

• if t is such that X^(x,v) = x* G <90 e , then S^x, t>) has the 
following jump discontinuity (in the velocity component only) 

(4.3) X{ +0 (x,v) = Xl (x,v) , V£ (^) = K(n x *)V t t (x,v) , 

where lZ{n x *) is the symmetry with respect to the hyperplane 
orthogonal to the inner unit normal n x * at x* G dQ e : in other 
words 

(4.4) K(n x .)S = Z-2(Z.n x .)n x ., ?GR D . 

Notice that this dynamics is mechanically reversible, i.e. for all t G 
R, one has 

S e _ t (x, v) = S e t {x, -v) , (x, v)en e x S - 1 . 

For each t G R, the map S% is invariant under Z D -translations, in the 
following sense: for all k G Z D and (x,v) G f2 e x S D_1 , one has 

5 t e (x + k,v) — (Xt(x, v) + fc, V t e (x, v)) = St(x, v) + (k, 0) . 
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Hence S% defines a one-parameter group — still denoted by — on 
the quotient space T e x S D_1 , where T e = Q e /Z D , and the restriction 
of the measure dxdv (defined on T D x S D_1 ) to T e x S D_1 is invariant 
under Sf. 

Starting from the billiard flow S% on T e x S ^ 1 , we define a unitary 
group §1 on L 2 (T t x S D_1 , dxdv) by the formula 

(4.5) Slf(x, v) = f(SU(x, v)) , (x, v) G T e x S D ~ l , t G R . 

Next we turn to the case of the same billiard table, but with ab- 
sorbing obstacles. By this, we mean that whenever a particle hits the 
boundary of Q e , it disappears. Equivalently, one may think of Q t as a 
sieve, with particles falling into the holes (the components of fl^). 

The analogue of in the case of absorbing obstacles is the contrac- 
tion semigroup T t e defined on L 2 (Y t x S D_1 , dxdv) by the formula 
(4.6) 

T^g(x, v) = g(x - tv, v)l Tr{x / e ^ v ) >t / e , (x, v ) G T t x S D 1 , t G R , 

where r r is the free path length defined in (12. ip — recall that r^e 1 ^ ^ 1 ^. 

Finally, L 2 (T t x S D ^ 1 ,dxdv) can be embedded as a subspace of 
L 2 (T D x S D_1 , dxdv ), identifying each function defined a.e. on T e x 
S D_1 with its extension by in the complement of T e x S ^ 1 . 

The maps and are related by the following elementary inequal- 
ity: for each / G L 2 (Y e x S ^ 1 , dxdv) and each t > 0, 

(4.7) / > a.e. implies that S%f(x,v) > Tjrf(x,v) a.e. 

on T e x S D_1 . Indeed, the formulas (IP]) . (jOj) . and imply that 

T tf( x > v ) = f( x -tv,v)l Tr{x/ei . v)>t/€ = f(S e _ t (x,v))l Tr{x/ei _ v)>t/€ 

<f(St t (x,v)) = Stf(x,v). 

Finally, we give the PDE interpretation of the operators and T t e . 
The function F e (t,x,v) = S^f(x,v) is the solution of 

d t F e + v ■ V x F e = , (x, v) G T e x S - 1 , 

F e (t, x, f ) = F 6 (t, x, TZ{n x )v) , (x, u) G <9T e x S D_1 , 
F\ = f 

(where TZ(n x ) designates the reflection with respect to the hyperplane 
orthogonal to n x : see fla. (14.41) ). while the function G t (t,x,v) = 
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T^g(x,v) is the solution of 

d t G e + v • V x G e = , (x,«)6T £ xS w , 
Gr e (£, x, v ) = , x E dT e , ?; • n x > , 

G 4=o = 9- 

4.2. Main result. The paper by Lorentz [15] described in the intro- 
duction suggests the following question, in the case of space dimension 
D = 3: 

"Let f in e L 2 (T 3 x S 2 ). Does S e t (f in l TtXS 2) (or any subsequence 
thereof) converge in L°°((0, +oo) x T 3 x S 2 ) weak-* as e — > to the 
solution / of fll.ip on T 3 x S 2 with N at = 1, r at = r#, and with initial 

data f n ?" 

The answer to that question is negative, as shown by the following 
theorem. 

Theorem 4.1. Assume that the space dimension satisfies D > 2. 
There exist initial data f m G L 2 (T D x S D_1 ) such that, for any a > 
and any function k 6 C(S D ~ 1 x S D_1 ) satisfying A3.3\) , no subsequence 
of Sl(f in l TeXS D-i) converges in L°°((0,+oo) x T D x S^" 1 ) weak-* to 
the solution of 

(4.8) (d t +v-V x )f{t,x,v) = a [ k{v,w)(f(t,x,w)- f(t,x,v))dw 

on T D x S -0-1 , rai/i initial data 

(4-9) /L =0 = / m - 

4.3. Proof of Theorem 14.11 We define 

$ e (t,2C,u) = l e r r (x/e,-v)>t ■ 

Since < $ e < 1, the sequence $ e is relatively weakly-* compact in 
L°°(R + x T D x S D_1 ). From now on, we consider a subsequence of 
e = 1/n, denoted e', so that $ e / -> $ in L°°(R+ x T D x S^" 1 ) weak-*, 
and denote r' = r*e /ly ' £>-1 . Then, the function $ is independent of x, 
i.e. 

(4.10) $ = $(t,v), 

as shown by the following classical lemma. 

Lemma 4.2. Let Z be a separable locally compact metric space, en- 
dowed with a Borel measure m. Let u n (x, z) = U n (nx, z), where (Z7 n ) n >i 
zs a bounded sequence of elements of L°°(T D x Z). Any weak-* limit 
point of the sequence (w n )n>i os n — > +oo zs a function of z alone (i.e. 
independent of x). 
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We postpone the proof of Lemma 14.21 until the end of this section. 
In the sequel, we consider an arbitrary initial data p = p(x) G 
L^iT ) that is independent of v and a.e. nonnegative. We define 

(4.11) f e (t,x,v) = 4 e (pl TeXSD -x) = p(Xl t {x,v))lr e (x) 
and 

(4.12) g e {t,x,v) = T t e (pl TeXSD - 1 ) =p(x-tv)$ e (t,x,v)lr t (x). 

Since lx e — > 1 a.e. on T D x S ^ 1 and |1tJ < 1, one has, by dominated 
convergence, 

(4.13) g e > -> p(x - tv)$(t, v) in L°°(R+ x T D x S ' 1 ) weak-* . 

On the other hand, (14.111) shows that ||/ £ ||z°° < ||p||l°°; therefore, possi- 
bly after extraction of a subsequence (still denoted e') , one has f e i — > f 
in L°°(R+ x T D x S 1 ^ 1 ) weak-*. We recall that p > a.e. on T D . 
Therefore, because of the inequality (14 .7p and of the weak-* limit f)4. 13j) . 
one has 

(4.14) f(t,x,v)>p(x-tv)$(t,v). 

In particular, 

(4.15) r 

f{t,x,v) 2 dxdv> II p(x — tv) 2 Q(t,v) 2 dxdv 

p{yfdy [ $(t,v) 2 dv 



> \\p\\h U i Ht,v)dv^j 



2 



where the last inequality follows from Jensen's inequality. By assump- 
tion, $ e / -> $ in L°°(R+ x T D x S^" 1 ) weak-*; because of (OD|l . one 
has 

$ £ /(t, x,v)dxdv —> / f)(if in L°°(R + ) weak-*. 



T D xS D-l 

D-1 



Furthermore, by Theorem 12.11 

<& e '(t, x, v)dxdv 

= meas({(x,v) G T e / x S " 1 | v(x/e', -v) > t/e'}) 
= meas({(y,v) G F r / x S 13-1 | v(?/, -t>) > t/e'}) 

> t^d-t = rS 1 ! > for a11 * > V^f _1 • 
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Hence, the inequality (14.151) becomes 
(4.16) 

f(t,x,v) 2 dxdv ) > — ^zrlMU 2 , for alH > 1/rf _1 . 



T D xS D-l 



tri 



Assume that / is the solution to (|£B|l - (|£gj) with f in = p. By Theo- 



rem [3111 one has 



p(x)dx 



< ce 



-7*1 



PlU 2 



L^TDxSfl- 1 ) 

for all t > 0. In particular, for all £ > 0, one has 

1/2 

\2 



(4.17) 



) < / p(x)dx + ce ^WpWl 2 ■ 



f(t, x, v) dxdv 

T D xS D-l 

In conclusion, if / is the solution to (I4.8p - (I4.9I) . then the initial data 
p (assuming it is not a.e. 0) must satisfy the inequality 



(4.18) 



< 



IpIU 2 



+ ce 



for alH > 1 /rl 



D-l 



At this point, we recall that this inequality holds for any arbitrary 
p G L°°(T D ) such that p > a.e.: hence the ratio ||p||z,i/||p||z,2 can be 
made arbitrarily small. For instance, one can choose p as follows: pick 
b, a bump function on R D satisfying 



1 llD 



0<b<l, supp(6)c[-i i 



For m G N*, define p to be the unique Z^-periodic function such that 



p i i : x h-> b(mi) . 

^ 2'2^ 



Then 



L 1 (T D ) 



m 



-D\ 



|l!(r d ) , 



i 2( T D) = m 



so that 



m 



~ D/2 -> as m 



-D/2i 



-oo . 



|l 2 (r d ) , 



|L 2 (T D ) 

Since the ratio HpIIlVIIpIU 2 can made arbitrarily small, the inequal- 
ity (14.180 would entail 



D-l 



< ce 



-7* 



for all i > 1 /r* 



D-l 



which is manifestly wrong for t large enough - - specifically, larger 
than the unique zero of the function t h- > Cie 7 * — cr®~H, denoted 
by t,(Ci,c,7,r»). 
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Hence, the assumption that, for each nonnegative p = p(x) G L°°(T D ), 
there exists a subsequence of (pl TeXS D-i) converging in L°° weak-* 
to the solution of (I4.8]l - (j4.9[) with f tn = p is wrong. This concludes the 
proof of Theorem 14.11 once Lemma 14.21 is proved. 

4.4. Proof of Lemma 14.21 Assume that 

u nq — > u in L°°(T D x Z) weak-* as n q — > +oo . 
For each k G Z D , one has 

u nq (k,z)= / e~* 2nk ' x u nk (x,z)dx = U nq (k,z) ii n q \k , 

= otherwise. 

Since the assumed convergence entails 

u nq (k, ■) — > u(k, •) in L°°(Z) weak-* for each k G Z D , 
this shows that 

•) = unless k = . 
Hence m is independent of announced. 

5. Final remarks 

We have demonstrated the impossibility of representing the Boltz- 
mann-Grad limit of the periodic Lorentz gas by a linear Boltzmann 
equation. 

Notice that this impossibility results solely from the lower bound on 
the distribution of free path lengths in Theorem 12.11 This is in fact 
not too surprising since the probabilistic representation of equation 
( 11.11) involves in particular exponentially distributed jump times for 
the velocity process. 

Notice also the choice of initial data in Theorem l4.ll Since one knows 
that particles moving in appropriately chosen rational directions may 
not encounter any obstacles on the billiard table Z r , it may seem some- 
what surprising that the contradiction in Theorem 14.11 is obtained by 
considering isotropic initial data (i.e. initial data that are indepen- 
dent of the angle variable) instead of pencils of particles concentrated 
in phase-space on those rational directions that avoid all obstacles. 
In fact, the contribution of such rational directions is already taken 
into account in the lower bound in Theorem 12.11 Besides, the case 
of isotropic initial densities is somewhat more natural in the context 
of Theorem 14.11 as it corresponds to local equilibria for the Lorentz 
kinetic model (11. ip . 
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As for the hydro dynamic (diffusion) limit of the periodic Lorentz gas 
with finite horizon, it was proved in pQ that the case of an isotropic 
reflection law at the surface of each obstacle can be treated by PDE 
techniques, avoiding the heavy machinery from ergodic theory required 
to handle the case of specular reflection and developed by Bunimovich- 
Sinai — see [7] , and also [6] . 

In the case of the Boltzmann-Grad limit however, changing the col- 
lision process does not affect the result in Theorem 14.11 since the ob- 
struction to using the linear Boltzmann equation comes from particles 
that travel too far before encountering an obstacle for the first time. 
For such particles, the nature of the collision process is obviously of 
no importance. Hence Theorem 14.11 holds verbatim if one replaces 
St(f m l TtXS D-i) with the solution of 

d t f £ + vVJ e = 0, xGT e , veS - 1 , 
/ e (t)|^ = /C(/ e (t)| EL ), 

/ I fin-i 

h\ t= Q — J J-T.xS 15 - 1 > 

where 

£± = {(x, v) G dT e x S D_1 | ±vn x >0} 

and /C is any linear operator from L°°(S1) to L°°(£+) that preserves 
the cone of positive functions. 

By the same token, the same result as in Theorem 14. II holds without 
change if the obstacles are not assumed to be spherical, or even identi- 
cal, but instead of arbitrary shapes, provided that they can be included 
in balls of radius r = r*e D//( - D_1 ' ) centered at the points of eZ D . 

Acknowledgements. I express my gratitude to Profs. C. Bardos, C. 
Boldrighini and H.S. Dumas for helpful comments during the prepara- 
tion of this paper. 
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